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Abstract 

We  develop  a  local  refinement  algorithm  for  analysis-suitable  T-splines  which  does  not  produce 
excessive  propagation  of  control  points.  We  then  demonstrate  its  use  as  an  adaptive  framework 
for  isogeometric  analysis.  Analysis-suitable  T-splines  are  a  class  of  T-splines  which  are  linearly 
independent  and  form  a  partition  of  unity.  These  properties,  coupled  with  local  refinement,  make 
this  class  of  T-splines  appealing  as  a  basis  for  isogeometric  analysis. 
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1.  Introduction 

Isogeometric  analysis  has  emerged  as  an  important  alternative  to  traditional  engineering  de¬ 
sign  and  analysis  methodologies.  Isogeometric  analysis  was  introduced  in  [1]  and  later  de¬ 
scribed  in  detail  in  [2],  In  isogeometric  analysis,  the  smooth  geometric  basis  is  used  as  the 
basis  for  analysis.  Most  of  the  early  developments  in  isogeometric  analysis  focused  on  establish¬ 
ing  the  behavior  of  the  smooth  NURBS  basis  in  analysis.  It  was  demonstrated  that  smooth¬ 
ness  offers  important  computational  advantages  over  standard  finite  elements  [3,  4].  Areas 
of  application  of  NURBS-based  isogeometric  analysis  include  turbulence  [5,  6,  7,  8],  fluid- 
structure  interaction  [9,  10,  11,  12],  incompressibility  [13,  14,  15],  structural  analysis  [16,  3], 
plates  and  shells  [17,  18,  19,  20,  21],  phase-field  analysis  [22,  23],  large  deformation  with  mesh 
distortion  [24],  shape  optimization  [25,  26,  27,  28],  and  electromagnetics  [29],  This  success 
has  in  turn  stimulated  efforts  within  the  Computer  Aided  Geometric  Design  (CAGD)  commu¬ 
nity  to  develop  and  integrate  analysis-suitable  geometric  technologies  and  isogeometric  analy¬ 
sis  [30,  31,  32,  33,  34,  35,  36,  37], 

While  smoothness  is  an  important  consideration,  NURBS  are  severely  limited  by  their  tensor 
product  construction.  In  traditional  NURBS-based  design,  modeling  a  complicated  engineering 
design  often  requires  hundreds,  if  not  thousands,  of  tensor  product  NURBS  patches  which  are 
usually  discontinuous  across  patch  boundaries.  Also,  almost  all  NURBS  models  use  trimming 
curves.  For  these  reasons,  a  global  geometric  discretization,  based  on  NURBS,  is  usually  not 
suitable  as  a  basis  for  analysis. 
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Figure  1:  The  domain  f2  C  M2  for  a  bivariate,  cubic  ( p  =  3)  T-spline. 


T-splines  were  introduced  as  a  superior  alternative  to  NURBS  [38].  T-splines  can  model 
complicated  designs  as  a  single,  watertight  geometry.  Additionally,  NURBS  are  T-splines  so 
existing  technology  based  on  NURBS  extends  to  T-splines.  Any  trimmed  NURBS  model  can  be 
represented  by  a  watertight  trimless  T-spline  [39]  and  multiple  NURBS  patches  can  be  merged 
into  a  single  watertight  T-spline  [38,  40],  Unlike  NURBS,  T-splines  can  be  locally  refined  [41]. 
These  properties  make  T-splines  an  ideal  technology  for  isogeometric  discretizations. 

Initial  investigations  using  T-splines  as  a  basis  for  isogeometric  analysis  demonstrated  that 
the  T-spline  basis  possesses  similar  convergence  properties  to  NURBS  with  far  fewer  degrees- 
of-freedom  [42,  43].  Additionally,  T-splines  possess  a  natural  finite  element  structure  which  can 
be  integrated  seamlessly  into  existing  finite  element  frameworks  via  Bezier  extraction  [44,  45]. 
T-splines  have  since  been  applied  to  problems  in  fracture  and  damage  [46,  47],  and  shells  [19]. 

Analysis-suitable  T-splines  were  introduced  in  [48].  Analysis-suitable  T-splines  are  a  mildly 
restricted  subset  of  T-splines.  Analysis-suitable  T-splines  are  linearly  independent  [48]  and,  if  a 
minor  boundary  condition  constraint  is  honored,  form  a  partition  of  unity  [49].  In  this  paper,  we 
develop  a  highly  localized  refinement  algorithm  for  analysis-suitable  T-splines  which  meets  the 
demands  of  both  design  and  analysis. 

This  paper  is  organized  as  follows.  Basic  T-spline  concepts  are  reviewed  in  Section  2. 
Analysis-suitable  T-splines  are  then  described  in  Section  3.  Section  4  presents  a  local  refinement 
algorithm  for  analysis-suitable  T-splines.  The  behavior  and  effectiveness  of  this  algorithm  is  then 
demonstrated  in  Section  5.  This  paper  focuses  on  bicubic  T-spline  surfaces,  although  the  con¬ 
cepts  generalize  to  arbitrary  odd  degree.  T-splines  of  arbitrary  degree  are  discussed  in  [42,  50]. 

2.  T-spline  fundamentals 

We  present  a  brief  overview  of  fundamental  T-spline  concepts  focusing  on  those  ideas  re¬ 
quired  to  understand  local  refinement.  A  more  detailed  description  of  T-splines  from  an  isogeo¬ 
metric  analysis  perspective  is  presented  in  [45].  We  base  our  developments  on  the  the  physical 
domain  Sic®2  shown  in  Figure  1 .  Throughout  this  paper  we  use  ds  to  indicate  the  number  of 
spatial  dimensions.  In  all  cases,  the  polynomial  degree  p  is  3. 
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2.1.  The  T-mesh 

The  fundamental  object  of  interest  underlying  T-spline  technology  is  the  T-mesh,  denoted 
by  T.  For  surfaces,  the  T-mesh  is  a  mesh  of  polygonal  elements.  Each  polygonal  element  is 
either  a  quadrilateral  or  an  element  with  quadrilateral  shape  where  one  or  more  edges  is  split  by 
T-junctions.  A  T-junction  is  analogous  to  a  “hanging  node”  in  finite  elements.  A  control  point, 
P a  €  W1k  ,  and  control  weight,  wa  €  R,  where  the  index  A  denotes  a  global  control  point 
number,  is  assigned  to  every  vertex  in  the  T-mesh.  A  mesh  for  the  domain  O  in  Figure  1  is  shown 
in  Figure  2a.  The  black  and  red  circles  are  T-mesh  vertices,  or,  equivalently,  control  points.  The 
T-junctions  in  Figure  2a  are  the  red  circles  P16,  Pn,  P28,  and  P2g. 

To  define  a  basis,  a  valid  knot  interval  configuration  must  be  assigned  to  the  T-mesh.  A 
knot  interval  [51]  is  a  non-negative  real  number  assigned  to  an  edge.  A  valid  knot  interval 
configuration  requires  that  the  knot  intervals  on  opposite  sides  of  every  element  sum  to  the  same 
value.  A  valid  knot  interval  configuration  for  the  T-mesh  in  Figure  2a  is  shown  in  Figure  2b. 


(a) 


(b) 


Figure  2:  The  T-mesh  and  knot  interval  configuration  defining  the  bicubic  T-spline  geometry  in  Figure  1.  (a)  The  T- 
mesh  defining  the  bicubic  T-spline  geometry  in  Figure  1.  The  red  circles  (16, 17,  28,  29)  are  T-junctions.  The  indexing 
identifies  the  T-mesh  control  points,  (b)  A  valid  knot  interval  configuration  for  the  bicubic  T-mesh  in  (a).  The  triangles 
correspond  to  a  knot  interval  of  0,  the  squares  correspond  to  a  knot  interval  of  1,  the  pentagons  correspond  to  a  knot 
interval  of  2,  and  the  stars  correspond  to  a  knot  interval  of  4.  Notice  that  the  knot  intervals  along  opposing  sides  of  each 
T-mesh  element  sum  to  the  same  value. 


2.2.  The  T-spline  basis 

Once  a  valid  knot  interval  configuration  has  been  assigned  to  a  T-mesh,  a  T-spline  basis  can 
be  constructed.  For  every  vertex  in  the  T-mesh,  a  T-spline  basis  function  is  constructed.  To 
illustrate,  we  construct  the  T-spline  basis  function  associated  with  Pi6-  We  note  that  a  class 
of  T-splines  where  the  T-spline  blending  functions  do,  in  fact,  constitute  a  basis  is  described  in 
Section  3  . 
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Figure  3:  Constructing  the  local  knot  interval  vectors  corresponding  to  T-mesh  vertex  Pi6- 


2.2.1.  Local  knot  interx’al  vectors 

A  T-spline  basis  function  is  constructed  from  knot  interval  sequences  inferred  from  the  T- 
mesh  in  the  neighborhood  of  the  associated  vertex.  These  knot  interval  sequences  are  called 
local  knot  interval  vectors.  A  local  knot  interval  vector  is  a  sequence  of  knot  intervals,  AS  = 
{A£l,  A£2,  •  ■  ■ ,  A£4},  such  that  A&  =  £i+i  - 

For  every  vertex.  A,  in  the  T-mesh,  we  construct  a  set  of  local  knot  interval  vectors,  AS4  = 
{AS^}?=1,  where  AS^  =  {A£^  1 ,  Af'A  2, . . . ,  A£a  5},  by  marching  through  the  T-mesh  in 
each  topological  direction,  starting  at  vertex  A.  until  2  vertices  or  perpendicular  edges  are  in¬ 
tersected.  At  each  intersection,  the  knot  interval  distance  traversed  since  the  last  intersection  is 
placed  in  the  local  knot  interval  vector.  If  a  T-mesh  boundary  is  crossed  before  2  knot  intervals 
are  intersected  it  is  common  to  set  the  remaining  knot  intervals  to  zero.  This  creates  an  open  knot 
vector  structure  along  the  boundary  of  the  T-mesh. 

A  set  of  local  knot  vectors,  H  .4  =  {S^}?=1,  where  ElA  =  l5  £A  2, . . . ,  £A  5},  can  be 

derived  from  A34  by  selecting  an  origin  O  £  R.  and  setting  £A  ■  =  O  +  YjL= 2  k- i-  Since 

a  knot  in  a  local  knot  vector  corresponds  to  T-mesh  topology  (vertices  or  edges)  encountered 
during  basis  function  inferrence  the  T-mesh  topology  and  knot  interval  configuration  determine 
the  knot  structure  of  the  underlying  T-spline  space.  In  the  context  of  T-splines,  the  term  knot  is 
often  also  used  to  refer  to  the  underlying  T-mesh  vertex  or  edge.  In  Figure  3,  the  knot  intervals 
used  to  construct  the  local  knot  interval  vectors  for  T-mesh  vertex  Pig  are  shown.  The  knot 
interval  vectors  for  P4g  are  given  by 


AS 


16  — 


1,1,  2,0 
0,2,  2,  2 


2.2.2.  The  local  basis  function  domain 

Using  AS^i  we  define  a  local  basis  function  domain,  LLa  C  R2,  as 

2 

^2 A  =  (^)  ^Ai 


i— 1 

4 


(1) 


where  il\  =  [0,  1  ^A,j]  C  K.  A  coordinate  system,  £A  =  (£^,  fA)  =  (£4,  pA),  called  the 

basis  coordinate  system ,  is  assigned  to  each  local  basis  function  domain. 

We  note  that  in  addition  to  the  basis  coordinate  systems  in  a  T-mesh,  it  is  often  desirable  to 
establish  larger  knot  coordinate  systems  for  a  subset  of  the  knot  structure  of  a  T-mesh.  Using  a 
knot  coordinate  system,  multiple  basis  functions  can  be  compared  in  a  common  coordinate  sys¬ 
tem.  Knot  coordinate  systems  are  used  when  computing  the  elements  of  the  refinement  operator, 
M,  as  described  in  Sections  2.4.2  and  4.3. 


2.2.3.  T-spline  basis  functions 

Over  each  local  basis  function  domain  f lA  we  define  a  single  T-spline  basis  function,  N A  : 
Qa  — >  R+  U  {0}.  This  is  done  by  forming  the  tensor  product  of  the  univariate  B-spline  basis 

functions  j NA3(^A  |  S^)|  as 

)  i— 1 

2 

iM£A|HA)  =  n^3(&lSA)-  (2) 

i- 1 


The  univariate  B-spline  basis  function,  NAP  :  Q.lA  — »•  M+  U  {0},  is  defined  using  a  recurrence 
relation,  starting  with  the  piecewise  constant  ( p  =  0)  basis  function 


(a  1  Ca,  1,  Ca,  2) 


1  ifa,i<a<a,2 

0  otherwise 


(3) 


where  tfA  k  is  the  kth  knot  value  in  the  local  knot  vector  d3A.  For  p  >  0,  the  basis  function  is 
defined  using  the  Cox-de  Boor  recursion  formula: 


N'7  7a  7a-7.2 . ^,,.2) 


pi  _  pi 

^A,l  7V7-2,p— 1 

pi  _  pi  VA 

^A,p+1  1 


{Pa  I  C 


A, If 


■  £>A,p+l 


) 


+ 


Ci  _  pi 

^ A .p-\-‘2  T.ji.p—1  (pi  I 

Ji  _  pi  lyA  WlM,2’ 

^yt,p+2  ^A,2 


A,p+2 


)• 


(4) 


2.3.  Bezier  element  construction 

A  Bezier  element  is  a  region  of  the  T-spline  surface  in  physical  space  bounded  by  knot  lines. 
Each  knot  line  in  physical  space  is  the  image  of  a  line  of  reduced  continuity  in  at  least  one  T- 
spline  basis  function.  We  call  the  collection  of  Bezier  elements  the  Bezier  mesh.  The  existence  of 
T-junctions  and  zero  knot  intervals  usually  results  in  there  not  being  a  one-to-one  correspondence 
between  T-mesh  elements  and  Bezier  elements. 

An  elemental  T-mesh,  Teiem,  is  formed  by  augmenting  T  with  all  images  of  basis  function 
knot  lines  that  do  not  already  correspond  to  an  edge  in  T,  then  eliminating  all  elements  for  which 
the  knot  interval  sum  on  any  side  is  zero.  The  elements  of  Teiem  are  in  one-to-one  correspon¬ 
dence  with  the  Bezier  elements.  Figure  4a  shows  Telem  for  the  T-mesh  in  Figure  2a.  Dashed 
lines  represent  the  edges  which  have  been  added. 

To  construct  appropriate  finite  element  paraphernalia  for  T-splines  we  use  T and  Bezier 
extraction  [45,  44],  Bezier  extraction  builds  a  linear  operator  for  each  Bezier  element.  The 
linear  transformation  is  defined  by  a  matrix  referred  to  as  the  extraction  operator.  The  extraction 
operator  maps  a  Bernstein  polynomial  basis  defined  on  Bezier  elements  to  the  global  T-spline 
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(a) 


(b) 


Figure  4:  The  elemental  T-mesh,  T eiern,  and  extracted  Bezier  mesh  in  physical  space  corresponding  to  the  T-mesh  in 
Figure  2a.  (a)  The  elemental  T-mesh,  T eiern.  The  dashed  lines  are  edges  which  have  been  added  to  T.  Only  the  non-zero 
parametric  area  elements  (the  shaded  elements)  are  included  in  T eiern.  The  indexing  identifies  the  elements  in  T eiern. 
(b)  The  extracted  Bezier  mesh  in  physical  space.  Each  element  in  (a)  corresponds  to  a  Bezier  element  in  physical  space. 
The  indexing  identifies  the  Bezier  elements. 


basis.  The  transpose  of  the  extraction  operator  maps  the  control  points  of  the  global  T-spline  to 
the  Bezier  control  points.  The  idea  is  illustrated  in  Figure  5  for  a  B-spline  curve.  This  element 
form  can  then  be  integrated  into  existing  finite  element  frameworks  in  a  straightforward  manner. 
See  [45]  for  additional  details. 

The  extracted  Bezier  mesh  in  physical  space  for  the  T-mesh  in  Figure  2a  is  shown  in  Fig¬ 
ure  4b.  Notice  the  one-to-one  correspondence  between  elements  in  Teiern  in  Figure  4a  and  the 
Bezier  elements  in  Figure  4b.  On  the  other  hand,  a  single  T-mesh  element  may  correspond  to 
multiple  Bezier  elements,  and  other  T-mesh  elements  may  not  correspond  to  any  Bezier  element. 
To  demonstrate,  the  Bezier  elements  corresponding  to  the  T-mesh  in  Figure  2  are  shown  in  Fig¬ 
ure  4b.  Notice  that  Bezier  elements  3  and  6  correspond  to  the  same  T-mesh  element,  as  do  Bezier 
elements  9  and  12. 

2.4.  T-spline  local  refinement 

The  set  of  all  T-splines  with  the  same  T-mesh  topology,  T,  and  knot  interval  configuration 
is  called  a  T-spline  space  [41].  We  denote  a  T-spline  space  by  T,  where  the  number  of  T-spline 
control  points  in  T  is  n.  While  the  notation  T1  C  T2  will  be  used  in  the  conventional  set- 
theoretic  sense,  the  notation  T1  C  T2  will  indicate  that  T2  can  be  created  by  adding  vertices  and 
edges  to  T1,  and  appropriately  modifying  the  knot  intervals  on  any  edges  which  are  split.  In  the 
context  of  finite  element  analysis,  vertices  and  edges  are  usually  added  by  subdividing  T-mesh 
elements. 
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Q  =  CTP 


Figure  5:  Bezier  extraction  for  a  B-spline  curve.  B-spline  basis  functions  and  control  points  are  denoted  by  N  and 
P,  respectively.  Bernstein  polynomials  and  control  points  are  denoted  by  B  and  Q,  respectively.  The  curve  T(£)  = 
PtN(£)  =  QrB(£).  The  extraction  operator  can  be  localized  to  the  individual  elements.  Note  that  the  Bernstein 
basis  is  the  same  for  each  element.  Formation  of  element  arrays  can  thus  be  standardized;  see  [45]  for  further  details. 


A  central  contribution  of  this  paper  (see  Section  4.2)  is  to  identify  conditions  under  which 

T1  c  T2  T1  c  r2.  (5) 

If  T1  C  T2,  T1  and  T2  are  said  to  be  nested  and  T1  is  a  local  refinement  of  T1  ■  In  Figure  6, 
a  T-mesh,  T1,  (solid  circles  and  lines)  and  corresponding  T-spline  space,  T1,  are  locally  refined 
through  the  addition  of  control  points  and  edges  (hollow  circles  and  dashed  edges).  In  this  case, 

T1  C  T2  ->  T1  C  T2. 

2.4.1.  Basis  function  refinement 

In  basis  function  refinement  [41],  knots  are  added  to  the  local  knot  vector  of  a  cubic  B-spline 
basis  function,  N(£  |  3),  where  3  =  {Ci,  C2>  •  •  • ,  Cs};  to  form  a  knot  vector,  3,  of  length  to.  N 
can  then  be  written  as  a  linear  combination  of  the  to  —  4  B-spline  basis  functions  defined  over 
substrings  of  length  5  in  3. 

For  the  case  to  =  6,  N(£  |  3)  is  split  by  inserting  a  single  knot  C  into  3  where  &  <  e  <  e  i+i- 
This  splits  the  basis  function  into  two  scaled  basis  functions: 

|  £l,  •••  ,£5)  =  aN(£  |  Cl;  •  •  •  )Ci:C;6+ 1;  •  •  •  >  60  +  bN(£  |  ^2,  ■  •  •  >  £»>  C;  Ci+1;  •  -  •  ■> &)  (6) 

where 

f  5  f°r  k  <  Ci 

a  =  {  ^  (7) 

1 1  for  k  >  Cr 
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Figure  6:  A  T-mesh,  T,  (solid  circles  and  lines)  and  T-spline  space,  T”1,  is  locally  refined  through  the  addition  of  control 
points  and  edges  (hollow  circles  and  dashed  edges).  In  this  case,  T1  C  T2  — »■  'T1  C  7~2. 


and 

f  |5~i  for  fc  >  £2 

6  =  ^  S  .  (8) 

\l  for  k  <  £2 

The  refinement  equations  for  the  case  m  >  6  can  be  derived  through  repeated  application  of 
these  equations. 

A  T-spline  basis  function,  N(£  \  S),  can  undergo  knot  insertion  in  either  parametric  direction 
by  inserting  a  knot  into  the  corresponding  local  knot  vector  and  then  applying  the  refinement 
equations.  This  results  in  two  scaled  T-spline  basis  functions  which  sum  to  the  original.  Further 
knot  insertion  into  these  resultant  scaled  basis  functions  yields  a  set  of  scaled  basis  functions 
which  also  sum  to  the  original. 

2.4.2.  The  refinement  operator  M 

if  r1  C  T2,  each  T-spline  basis  function,  N\  £  7-1,  can  be  expressed  uniquely  as  a  linear 
combination  of  the  T-spline  basis  functions,  1V|  £  T2,  as 

n2 

N\=J2  mA,BN%  (9) 

B=1 

where  n2  is  the  number  of  control  points  in  T2  and  the  toa,s  are  determined  by  knot  insertion 
as  described  in  Section  2.4.1.  This  relationship  can  be  written  in  matrix- vector  notation  as 

N1  =  MN2  (10) 

where  N1  =  {N±,  N%,  ■  ■  ■ ,  N^}7  is  the  column  vector  of  T-spline  basis  functions,  N\  £  T1, 
N2  =  {N’i,  JVf , . . . ,  N^}1  is  the  column  vector  of  T-spline  basis  functions,  Ng  £  T2,  and  M 
is  an  rii  x  n2  matrix  with  elements  rnA  [S.  We  call  M  the  refinement  operator. 
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3.  Analysis-suitable  T-splines 


Analysis-suitable  T-splines  form  a  practically  useful  subset  of  T-splines.  Analysis-suitable 
T-splines  maintain  the  important  mathematical  properties  of  the  NURBS  basis  while  providing 
an  efficient  and  highly  localized  refinement  capability.  All  T-splines  possess  the  following  prop¬ 
erties: 

•  The  basis  constitutes  a  partition  of  unity  [49]  (see  Section  3.4.) 

•  Each  basis  function  is  non-negative. 

•  An  affine  transformation  of  an  analysis-suitable  T-spline  is  obtained  by  applying  the  trans¬ 
formation  to  the  control  points.  We  refer  to  this  as  affine  covariance.  This  implies  that  all 
“patch  tests”  (see  [52])  are  satisfied  a  priori. 

•  They  obey  the  convex  hull  property. 

•  They  can  be  locally  refined. 

While  most  (but  not  all  [53])  T-splines  are  also  linearly  independent,  analysis-suitable  T- 
splines  are  always  linearly  independent  for  any  choice  of  knot  intervals  [48].  Analysis-suitable 
T-spline  spaces  are  defined  over  a  mildly  restricted  set  of  allowable  T-mesh  topologies.  This 
topological  restriction  can  be  described  elegantly  in  terms  of  T-junction  extensions. 

3.1.  T-junction  extensions 

A  T-junction  extension  is  normally  composed  of  a  face  and  edge  extension  (if  one  exists.) 
We  define  face  and  edge  extensions  to  be  closed,  directed  line  segments  which  originate  at  a 
T-junction.  An  extended  T-mesh  Ti  xt  is  formed  by  adding  all  T-junction  extensions  to  a  T-mesh 
T.  The  extended  T-mesh,  Text ,  for  the  T-mesh  in  Figure  7a  is  shown  in  Figures  7b  and  7c, 
respectively.  The  dotted  black  arrows  are  face  extensions  and  the  dashed  red  arrows  are  edge 
extensions.  The  T-junctions  are  denoted  by  large  circles. 

Figure  7b  shows  the  T-junction  extensions  in  the  index  space  [42]  of  the  T-mesh.  The  index 
space  is  created  by  plotting  the  knots  in  the  T-mesh  at  equally  spaced  intervals,  regardless  of 
their  actual  values.  The  index  space  point  of  view  is  useful  for  developing  algorithms,  as  well 
as  for  building  intuition.  For  example,  in  the  index  space  it  is  easy  to  identify  the  knot  lines  at 
which  the  support  of  any  given  function  will  begin  or  end.  Additionally,  the  direction  of  traversal 
and  orientation  of  a  T-junction  extension  can  be  uniquely  established  in  the  index  space  of  the 
T-mesh. 

Figure  7c  shows  the  set  of  T-junction  extensions  drawn  on  the  T-mesh  in  physical  space.  We 
often  drop  the  distinction  between  the  index  and  physical  space  representation  for  T-junction 
extensions  and  use  the  physical  space  representation. 

A  T-junction  extension  is  formed  in  a  manner  similar  to  what  was  described  for  the  construc¬ 
tion  of  local  knot  interval  vectors  in  Section  2.2.1.  A  face  extension  is  created  by  marching  from 
the  T-junction,  in  the  direction  of  a  missing  edge,  until  two  perpendicular  edges  or  vertices  are 
intersected.  The  direction  of  an  extension  is  always  away  from  its  T-junction.  An  edge  extension 
is  then  formed  only  if  an  edge  is  attached  to  the  T-junction  in  the  opposite  direction.  If  so,  the 
extension  is  formed  by  marching  in  the  opposite  direction  of  the  face  extension  until  the  edge’s 
opposite  vertex  is  encountered.  Since  T-junction  extensions  are  closed  line  segments,  a  horizon¬ 
tal  and  vertical  extension  can  intersect  either  on  the  interior  of  both  extensions  or  at  the  endpoint 
of  one  extension  or  both  extensions. 
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As  an  example,  consider  T-junction  E  in  Figure  7b.  The  face  extension  points  to  the  right 
and  intersects  the  two  vertical  edges  corresponding  to  indices  5  and  6  along  the  bottom.  Since 
T-junction  E  is  connected  to  an  edge  in  the  direction  opposite  the  face  extension,  we  also  form  an 
edge  extension  along  that  edge.  The  edge  extension  points  to  the  left  and  intersects  the  vertical 
edge  corresponding  to  index  2  along  the  bottom.  The  T-junction  extension  for  T-junction  E  is 
composed  of  the  face  and  edge  extension. 

3.2.  The  extension  graph 

Intersecting  T-junction  extensions  in  an  extended  T-mesh  T ext  can  be  visualized  using  an 
undirected  graph.  We  call  this  graph  the  extension  graph  and  denote  it  by  E(Text).  Each  node  in 
E  corresponds  to  a  single  T-junction  extension  in  Text.  If  two  extensions  in  Text  intersect  then 
an  edge  is  drawn  between  the  corresponding  nodes  in  E.  The  extension  graph  for  the  extended 
T-mesh  in  Figure  7b  is  shown  in  Figure  8a.  In  this  case  there  are  five  intersections  represented 
by  the  five  edges  in  the  graph. 

3.3.  Analysis-suitable  definition 

An  analysis-suitable  T-spline  is  one  whose  extended  T-mesh  is  analysis-suitable.  An  analysis- 
suitable  extended  T-mesh  is  one  where  no  T-junction  extensions  intersect.  In  other  words, 
E(T f  xt  )  is  an  empty  graph  (no  edges  in  the  graph).  We  denote  an  analysis-suitable  T-spline  space 
by  Ts  and  analysis-suitable  T-mesh  by  Ts.  The  T-mesh  in  Figure  7a  is  not  analysis-suitable.  This 
can  be  seen  by  inspecting  the  extension  graph  in  Figure  8a  which  has  five  edges.  By  adding  the 
dashed  edges  in  Figure  8b  the  T-mesh  becomes  analysis-suitable  because  its  extension  graph  is 
empty. 


3.4.  Partition  of  unity 

An  analysis-suitable  T-spline  basis  forms  a  partition  of  unity.  In  other  words,  X]  4=1  ^ a  =  1. 
The  partition  of  unity  property  is  important  for  both  geometry  and  analysis  because  it  assures 
affine  covariance  and  exact  satisfaction  of  all  patch  tests.  In  the  context  of  T-splines,  the  partition 
of  unity  property  can  be  described  in  two  ways.  The  equation  for  any  T-spline  surface  is 


J2a=i  P awaNa 

ELi  waN A 


(11) 


where  the  P  |  are  control  points,  w  a  are  weights,  and  Na  are  blending  functions.  This  equation 
can  also  be  written  as 


t = e tnn\  ^ 


(12) 


in  which  case  the  Ra  may  be  rational  T-spline  blending  functions.  It  is  clear  that  the  R a  always 
sum  to  one,  regardless  of  the  choice  of  wa-  In  general,  affine  covariance  only  requires  that 
the  Ra  form  a  partition  of  unity,  not  the  Na-  However,  it  is  shown  in  [49]  that  for  an  analysis- 
suitable  T-spline  with  all  wa  =  1,  Ra  =  Na-  Thus,  when  we  say  that  analysis-suitable  T-splines 
form  a  partition  of  unity,  we  mean  that  both  the  Na  and  Ra  sum  to  one.  This  stronger  notion  of 
partition  of  unity  is  also  a  property  of  NURBS. 
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Figure  7:  The  extended  T-mesh,  Text,  in  the  index  and  physical  space.  The  dotted  arrows  (in  black)  represent  face 
extensions  and  the  dashed  arrows  (in  red)  represent  edge  extensions.  The  T-junctions  are  denoted  by  large  circles,  (a)  A 
T-mesh  T  with  six  T-junctions,  (b)  The  extended  T-mesh  formed  from  T  and  the  T-junction  extensions  in  index  space. 
The  indexing  of  the  knot  structure  of  T  is  along  the  bottom  and  left,  (c)  The  extended  T-mesh  in  physical  space.  Note 
that  the  directionality  of  each  extension  is  always  determined  in  the  index  space  of  the  T-mesh. 
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(a) 


(b) 


Figure  8:  The  extension  graph,  E(Tea;£)-  (a)  The  extension  graph  E(T ext)  corresponding  to  T ext  in  Figures  7b  and  7c. 
The  five  edges  in  the  graph  correspond  to  the  five  intersections  between  T-junction  extensions  in  Text-  (b)  The  T-mesh 
in  Figure  7a  can  be  made  analysis-suitable  by  adding  the  bold  dashed  edges.  The  extension  graph  for  this  new  T-mesh  is 
empty  (no  edges.) 
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3.5.  The  analysis-suitable  elemental  T-mesh 

For  analysis-suitable  T-splines,  the  elemental  T-mesh,  T eiem  (see  Section  2.3),  can  be  formed 
by  simply  adding  the  face  extensions  to  T  [49],  This  greatly  simplifies  the  construction  of  T eiem 
since  it  is  not  necessary  to  inspect  the  knot  lines  in  the  T-spline  basis. 

4.  Local  refinement  of  analysis-suitable  T-splines 

The  T-spline  local  refinement  algorithm  presented  in  [41]  may  add  many  superfluous  control 
points  to  a  T-mesh  during  refinement.  Additionally,  using  this  algorithm  to  locally  refine  an 
analysis-suitable  T-spline  often  results  in  a  refined  T-spline  which  is  not  analysis-suitable. 

This  behavior  can  be  attributed  to  the  generality  of  the  algorithm,  which  is  designed  to  op¬ 
erate  on  any  T-spline.  No  assumptions  are  made  about  the  topological  characteristics  of  the 
underlying  T-mesh  or  space.  By  restricting  ourselves  to  analysis-suitable  T-splines,  however,  we 
can  leverage  the  structure  of  the  T-mesh  to  develop  a  simple  local  refinement  algorithm  which 
only  introduces  a  minimal  number  of  superfluous  control  points  and  preserves  the  properties  of 
an  analysis-suitable  space. 

4.1.  Analysis-suitable  nesting  theory 

We  say  T2ext  is  a  refinement  of  T \xt  (denoted  T lxtQT2xt)  if  T1  C  T2  and  no  face  extension 
endpoint  in  T f  xt  corresponds  to  a  point  in  the  interior  of  a  face  extension  in  T \xt  of  the  same 
directionality.  Then,  if  T}  and  T2  are  analysis-suitable  T-spline  spaces,  TlxtQT 2xt  — >■  T}  C 

[54],  Note  that  we  compare  face  extensions  in  the  index  space  of  T  jxt. 

Nestedness  between  two  T-spline  spaces,  T}  and  T2  (not  necessarily  analysis-suitable),  can 
be  visualized  using  a  simple  modification  of  the  extension  graph  described  in  Section  3.2.  We 
call  this  graph  the  coupled  extension  graph ,  E(T^~).  T^2,  called  a  coupled  extended  T-mesh, 
is  constructed  by  adding  the  face  extensions  of  T }.xt  to  T2ext. 

To  construct  i?(T^2)  we  augment  E{T2ext)  by  adding  an  additional  edge  to  the  graph  if  a 
T-junction  face  extension  endpoint  in  T2xf  is  in  the  interior  of  a  face  extension  from  T \xt  of  the 
same  directionality.  In  that  case,  a  “loop  edge”  is  drawn  from  the  corresponding  node  in  E(T2xt) 
to  itself  creating  a  small  loop  in  the  graph.  If  T2ext  is  analysis-suitable  (E(  T2xt)  is  empty)  then 
the  only  edges  which  exist  in  E( T^2)  are  “loop  edges.”  The  non-existence  of  loop  edges  is  a 
necessary  condition  that  TgXtCT2xi.  If  ^(T^2)  is  an  empty  graph  then  T}  C  T 2.  The  weight 
of  a  coupled  extension  graph  is  the  number  of  edges  in  the  graph  and  is  denoted  by  W{E).  The 
weight  of  a  node  is  the  number  of  edges  touching  the  node. 

Figure  9  illustrates  the  construction  of  a  coupled  extension  graph  and  its  use  in  determining 
nestedness.  An  analysis-suitable  extended  T-mesh  T \xt  is  shown  in  Figure  9a  and  T2xt  is  shown 
Figure  9b.  We  construct  the  coupled  extended  T-mesh  T^2  in  Figure  9c.  Notice  that  in  this  case 
the  only  face  extension  from  T \xt  which  is  visible  (a  light  gray  dotted  arrow  in  a  box)  corresponds 
to  T-junction  extension  B  in  Figure  9a.  This  indicates  that  the  endpoint  of  T-junction  extension 
C  in  T2X(  (see  Figure  9b)  is  in  the  interior  of  T-junction  extension  B  in  T\xt  (see  Figure  9a.) 
Figure  9d  shows  E{ T^2).  Since  the  graph  is  not  empty,  T}  ^  T2.  In  fact,  the  edges  between 
different  nodes  means  the  underlying  extension  graph  E(T'2xt)  is  not  empty  which  implies  that 
T2  is  not  analysis-suitable.  The  loop  edge  which  begins  and  ends  at  node  C  indicates  that 
T1  <iT2 

ext  y—  ext’ 
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(C)  (d) 


Figure  9:  Determining  nestedness,  (a)  An  analysis-suitable  extended  T-mesh  Tj  t.  (b)  An  extended  T-mesh  T2  +  (not 
analysis-suitable),  (c)  Superimposing  the  face  extensions  from  (a)  on  T2xt  to  form  T^2.  Notice  that  the  only  face 
extension  from  T \xt  which  is  visible  is  in  the  light  gray  box.  (d)  The  coupled  extension  graph  for  (c).  The  graph  contains 
edges  which  implies  that  T g1  £  7~2 . 
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4.2.  A  local  refinement  algorithm 

The  example  in  Figure  9  motivates  a  simple  approach  to  local  refinement.  First,  create  T2  D 
Tj .  As  in  standard  finite  element  analysis,  this  is  usually  done  by  subdividing  a  set  of  T-mesh 
elements  in  T2.  The  T-mesh  elements  are  often  selected  so  as  to  reduce  error  in  the  finite  element 
solution.  If  E(Tl^2)  is  not  empty,  we  must  add  some  additional  control  points  and  edges  to  T2 
to  cause  E{T\^t2)  to  be  empty. 

To  illustrate.  Figure  10a  shows  the  T-mesh  T2  from  Figure  9b,  where  T}  £  T2  eventhough 
Tj  C  T2.  This  can  be  seen  by  inspecting  the  coupled  extension  graph  in  Figure  9d.  To  ensure 
nestedness,  additional  refinement  of  T2  must  be  performed.  Figures  10b  through  lOd  show 
three  possible  analysis-suitable  refinements  of  T2  where  the  resulting  coupled  extension  graph 
EiJle7t2)  is  empty.  The  dashed  edges  and  open  circles  in  Figure  10  are  T-mesh  edges  and 
vertices,  respectively,  added  during  local  refinement.  The  minimum  of  the  three  refinements  is 
shown  in  Figure  lOd,  where  only  6  vertices  and  8  edges  have  been  added.  We  note  that  Figure  10b 
is  a  NURBS  refinement. 

This  example  raises  the  question  of  how  to  devise  an  algorithm  for  automatically  finding 
the  fewest  additional  control  points  and  edges  that  will  cause  E{J]fi^t2)  to  be  empty.  Finding 
the  minimal  number  is  an  NP-hard  problem,  so  for  efficiency,  our  algorithm  uses  the  following 
greedy  strategy  that  only  provides  an  approximate  minimum.  The  following  steps  constitute  our 
analysis-suitable  local  refinement  algorithm: 

1.  Create  T2  DTj. 

2.  Using  T*  and  T2,  form  E( T^2). 

3.  Of  all  possible  insertion  edges,  add  one  into  T2  for  which  the  weight  of  the  resulting 
E(7l7t2)  is  smallest.  An  insertion  edge  has  a  vertex  that  is  a  T-junction  in  T2  and  the 
corresponding  node  in  the  coupled  extension  graph  has  non-zero  weight  (see  Figure  11.) 

4.  Repeat  Step  3  until  the  weight  of  E(Tl^2)  is  zero. 

5.  Compute  the  refinement  operator  M,  if  desired.  See  Section  4.3. 

Only  one  insertion  edge  is  inserted  during  each  iteration  of  the  refinement  algorithm.  Also, 
there  are  cases  for  which  the  weight  will  stay  the  same  or  increase  after  the  addition  of  the  optimal 
insertion  edge.  It  should  be  noted  that  the  algorithm  will  always  terminate,  because  in  the  limit, 
if  all  T-junctions  are  extended  all  the  way  to  a  boundary  edge,  a  NURBS  is  created. 

We  first  demonstrate  Steps  1  through  4  of  analysis-suitable  local  refinement  on  a  simple 
example.  Step  5  is  explained  in  detail  in  Section  4.3.  We  begin  with  the  analysis-suitable  T-mesh 
T2  shown  in  Figure  12a.  In  Figure  12b,  T'  is  refined  by  subdividing  several  T-mesh  elements. 
The  coupled  extension  graph  E( T^2)  can  then  be  constructed  from  the  coupled  extended  T- 
mesh  Tgj^2  as  shown  in  Figures  12d  and  12c,  respectively.  Notice  that  the  weight  of  the  graph 
is  4  so  nesting  is  not  assured  between  T2  and  T2 ■ 
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(C) 


(d) 


Figure  10:  Several  analysis-suitable  local  refinements  for  the  example  in  Figure  9.  The  dashed  edges  and  open  circles 
are  T-mesh  edges  and  vertices,  respectively,  added  during  local  refinement,  (a)  The  T-mesh  T2  from  Figure  9b.  (b)  18 
vertices  and  19  edges  have  been  added,  (c)  7  vertices  and  8  edges  have  been  added,  (d)  6  vertices  and  8  edges  have  been 
added. 
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(a) 


(b) 


Figure  1 1 :  T-mesh  edges  which  may  be  inserted  during  Step  3  of  the  local  refinement  algorithm,  (a)  The  T-mesh  T2 .  (b) 
The  thick  black  dashed  lines  are  edges  which  can  be  inserted  into  T2 .  Notice  that  each  of  these  edges  has  a  vertex  which 
is  a  T-junction  in  T2  and  the  corresponding  node  in  the  coupled  extension  graph  has  non-zero  weight  (see  Figure  9d.) 
The  thick  black  dashed  lines  are  called  insertion  edges.  The  thin  red  dashed  lines  are  examples  of  T-mesh  edges  which 
are  not  inserted  to  form  a  refined  T-mesh  because  they  do  not  satisfy  the  requirements  of  an  insertion  edge. 
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(a) 


(b) 


(c) 


(d) 


Figure  12:  Initialization  of  a  simple  analysis-suitable  local  refinement  example,  (a)  The  initial  analysis-suitable  T-mesh 
before  refinement,  (b)  Four  T-mesh  elements  are  subdivided  to  form  T2.  (c)  The  coupled  extended  T-mesh  Tg^2 
corresponding  to  (a)  and  (b).  (d)  The  coupled  extension  graph  ^(Tg^2).  The  graph  is  not  empty  so  nesting  does  not 
necessarily  hold. 
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We  now  begin  to  add  insertion  edges  to  T2  until  is  empty.  The  result  of  the  first 

iteration  of  the  algorithm  is  shown  in  Figure  13.  The  input  coupled  extension  graph  is  shown 
in  Figure  13a.  The  corresponding  set  of  insertion  edges  is  shown  in  Figure  13b.  Notice  that 
insertion  edges  are  only  created  for  T-junctions  where  the  weight  of  the  corresponding  node  in 
E{J\7t2)  is  nonzero.  Each  insertion  edge  is  then  added  to  T2  as  a  T-mesh  edge  and  the  change 
in  graph  weight  AW  and  total  resulting  graph  weight  W  are  computed  as  shown  in  the  Table  of 
Figure  13.  The  shaded  cells  are  the  AlTs  which  must  be  computed  during  this  iteration.  Since 
this  is  the  first  iteration  all  must  be  computed.  The  insertion  edge  with  the  largest  AW  is  then 
selected  and  inserted  into  the  T-mesh.  For  this  iteration,  insertion  edge  K  (the  bold  dashed  line 
in  Figure  13b)  is  inserted  into  T2.  Notice  that  in  this  case  insertion  edge  E  could  also  have  been 
selected.  Both  have  a  AW  of  2. 


(a)  Input  graph  (weight  =  4) 


(b)  Insertion  edges 
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Figure  13:  The  first  local  refinement  iteration  for  the  example  in  Figure  12.  (a)  The  input  coupled  extension  graph.  This 
graph  has  a  weight  of  4.  (b)  The  corresponding  insertion  edges.  AW  and  W  for  each  insertion  edge  are  shown  in  the 
Table.  AW  measures  the  change  in  graph  weight  after  the  edge  is  inserted  into  I  IT  is  the  total  graph  weight  after 
the  edge  is  inserted  into  T2.  The  shaded  cells  indicate  AW  values  which  were  computed  during  this  iteration.  Insertion 
edge  K  minimizes  the  graph  weight  and  is  inserted  into  the  T-mesh  as  a  T-mesh  edge.  Notice  that  in  this  case  insertion 
edge  E  could  also  have  been  selected.  Both  have  a  AW  of  2. 


The  second  and  final  iteration  of  the  algorithm  is  shown  in  Figure  14.  The  coupled  extension 
graph  is  shown  in  Figure  14a.  Notice  the  presence  of  the  new  T-mesh  edge  in  T2  which  cor¬ 
responds  to  insertion  edge  K  from  the  previous  step.  The  current  insertion  edges  are  shown  in 
Figure  14b.  There  are  now  only  three  since  the  previous  iteration  decoupled  nodes  A,  J,  and  K. 
The  current  values  of  AW  and  W  are  shown  in  the  Table  of  Figure  14.  None  of  the  AW  cells  are 
shaded  which  indicates  that  the  values  are  saved  from  a  previous  iteration  and  not  computed.  In 
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(a)  Input  graph  (weight  =  2) 


(b)  Insertion  edges 
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Figure  14:  The  second  and  final  local  refinement  iteration  for  the  example  in  Figure  12.  (a)  The  input  coupled  extension 
graph.  This  graph  has  a  weight  of  2.  (b)  The  corresponding  insertion  edges.  AW  and  W  for  each  insertion  edge  are 
shown  in  the  Table.  AW  measures  the  change  in  graph  weight  after  the  edge  is  inserted  into  T2.  W  is  the  total  graph 
weight  after  the  edge  is  inserted  into  T2 .  Since  no  A  W  cells  are  shaded  all  values  are  saved  from  the  previous  refinement 
step  (see  Figure  13.)  Insertion  edge  E  drives  the  graph  weight  to  zero  and  is  inserted  into  the  T-mesh  as  a  T-mesh  edge  to 
complete  the  refinement  process. 


general,  only  a  small  number  of  AfVs  need  to  be  computed  during  each  step  of  local  refinement. 
The  insertion  edge  with  the  largest  AW  is  E  which,  when  inserted,  drives  the  total  weight  of 
the  graph  to  zero  thus  terminating  the  refinement  process.  The  final  refined  T-mesh  is  shown  in 
Figure  15.  The  new  edges  are  the  dashed  lines.  If  desired,  the  refinement  operator  M  can  be 
computed  as  described  in  Section  4.3. 

4.3.  Computing  the  refinement  operator  M 

We  now  describe  how  the  elements,  of  a  refinement  operator,  M,  are  computed.  We 

recall  that  Steps  1  through  4  of  the  analysis-suitable  local  refinement  algorithm  in  Section  4.2 
guarantee  that  T}  C  Tfi  and  the  existence  of  M.  We  note  that  all  knot  comparisons  between 
basis  functions  are  made  in  a  common  knot  coordinate  system  defined  in  the  index  space  of 
(see  Section  2.2.2.) 

Before  proceeding,  we  define  several  important  index  space  concepts.  For  a  basis  func¬ 
tion,  N\  g  Tsk,  we  can  construct  the  index  vectors ,  =  {F Afcli=i>  where  k  = 

{ta\  ,  t^’ 2,  ■  ■  • ,  5}  and  *•  is  the  index  of  in  the  index  space  of  Tfc.  Using  r^fc  we 
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Figure  15:  The  final  refined  T-mesh  for  the  example  in  Figure  12.  The  new  edges  and  vertices  are  the  dashed  lines  and 
open  circles,  respectively. 


can  then  define  the  index  domain  SlJ \  k  C  M2  as 

2 

i= 1 


03) 


Where  nA,k  =  [Ta.I-Ta’s]  C  R' 

We  first  initialize  all  entries  in  M  to  zero.  Then,  for  each  NA  £  Ts[  and  NB  £  T)2  such 
that  fig  2  C  f VA  2,  we  insert  £lB  into  EA  if  ^  2,  i  =  1,2 ,  j  =  1, 2, . . . ,  5.  This 

constructs  refined  local  knot  vectors,  3,4  =  {Z\}2i=1.  We  then  apply  basis  function  refinement 
(see  Section  2.4.1)  using  3 4  to  generate  the  linear  combination 

m 

na  =  Y.CA*Nk-  (14) 

fc=i 

If  there  exists  an  such  that  =  NB,  then  we  set  rriA,B  =  CA,k ■  Note  that  the  index  vectors 
and  domains  for  N\  are  consUucted  using  the  index  space  of  T 2.  This  is  possible  since  T'  C  T2. 

To  illustrate.  Figure  16a  shows  a  T-mesh,  T*,  where  additional  vertices  and  edges  (open 
circles  and  dashed  lines)  have  been  added  during  the  first  four  steps  of  the  local  refinement 
algorithm  in  Section  4.2  to  form  T2.  As  a  result  T}  C  T)2.  The  knot  intervals  used  in  this 
example  are  shown  next  to  the  corresponding  edges.  The  index  space  representation  is  shown  in 
Figure  16b.  For  simplicity  we  choose  a  global  knot  coordinate  system  with  an  origin  at  (1,  0)  in 
the  index  space  of  T2.  We  note  that,  in  practice,  the  knot  coordinate  system  can  be  defined  by 
the  subset  of  knots  in  T2  which  define  N\  and  NB. 

We  now  compute  where  the  basis  functions  N\  £  T}  and  NB  £  Tg  are  associated 

with  the  vertices  labeled  A  and  B  in  Figure  16.  In  this  case,  the  index  vectors,  with  respect  to  the 
index  space  of  T2,  are 

"  1,2, 3, 5, 6 
0,1, 2, 4, 5 
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(b) 


Figure  16:  Computing  the  refinement  coefficient,  rriA,B •  (a)  The  T-mesh,  T*,  consists  of  the  solid  circles  and  lines. 
The  topology  (control  points  and  edges)  added  during  the  topology  phase  of  analysis-suitable  local  refinement  (Steps  1 
through  4  in  Section  4.2)  are  denoted  by  the  open  circles  and  dashed  lines,  respectively.  The  topology  phase  ensures 
that  Tg  Q  7?.  (b)  The  index  space  of  T^.  The  index  domain  2  =  [1»  6]  0  [0,  5]  is  the  union  of  the  dark  and 
lightly  shaded  rectangles  and  the  index  domain  Qg  2  =  [2,  6]  0  [0, 4]  is  the  lightly  shaded  rectangle.  The  origin  of  the 
common  knot  coordinate  sys  tern  is  (1,  0)  in  the  index  space.  The  knot  intervals  for  the  common  knot  coordinate  system 
are  shown  in  (a)  and  (b). 


and 


r_e,2  = 


2, 3, 4, 5, 6 

0,1, 2, 3, 4 


06) 


The  index  domain,  2  =  [1, 6]  ®  [0,  5],  is  the  union  of  the  dark  and  lightly  shaded  rectangles 
in  Figure  16b  and  the  index  domain,  iVn  2  =  [2, 6]  8  [0, 4],  is  the  lightly  shaded  rectangle. 

In  this  coordinate  system,  the  local  knot  vectors  for  N\  and  Nj}  are 


0,1,  2, 4, 4 
0,0,  2, 4, 4 


(17) 


and 


1,2, 3, 4, 4 
0,0,  2, 4, 4 


08) 


Obviously,  VLtb  2  C  nA  2,  and  since  tb  3  = 
the  corresponding  knots  3  =  3  and  4 
knot  vectors 


4  and  tb24  =  3  are  not  in  2  and  2  we  insert 
=  4  into  E\  and  EA.  This  results  in  the  refined  local 

0,1,  2,  3, 4, 4 
0,0,  2, 4, 4, 4 
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(19) 


Applying  the  refinement  equations,  (6)  -  (8),  to  H  4,  in  each  univariate  direction,  and  taking  the 
tensor  product  of  the  resulting  coefficients  results  in 


where 


and 


N\  —  Ca,\N\  +  CA,  2N2  +  Ca,3^3  +  CA,iNi 


=  3W>  + 


0,1,  2, 3, 4 

0,0, 2, 4, 4 

1,2, 3, 4, 4 
0,0,  2, 4, 4 


Since  iV3  =  Ng,  we  have  that  niA,B  =  j. 


(20) 

(21) 


(22) 


(23) 


5.  Applying  analysis-suitable  local  refinement 

We  now  explore  the  application  and  behavior  of  analysis-suitable  T-splines  and  local  refine¬ 
ment.  In  this  example,  we  use  analysis-suitable  local  refinement  to  transform  an  initial  coarse 
T-spline  ship  hull  design  into  an  analysis-suitable  model.  The  analysis-suitable  model  can  then 
be  used  directly  in  isogeometric  analysis  by  way  of  Bezier  extraction  [45].  We  note  that  this 
same  approach  can  also  be  used  as  an  adaptive  finite  element  solution  strategy.  A  demanding 
application  of  adaptive  T-spline  local  refinement  in  the  context  of  a  phase-field  fracture  model  is 
described  in  [55]. 

In  Figure  17,  we  schematically  illustrate  the  process  used  to  perform  analysis-suitable  local 
refinement.  First,  a  set  of  Bezier  elements,  generated  using  Bezier  extraction,  is  flagged  by  the 
user  or  finite  element  solver.  Next,  the  Bezier  elements  are  used  to  identify  corresponding  T-mesh 
elements.  We  recall  that  several  Bezier  elements  may  correspond  to  a  single  T-mesh  element  as 
described  in  Section  2.3  and  shown  in  Figure  17  on  the  left.  The  selected  T-mesh  elements  are 
then  refined  to  generate  T2  as  described  in  Section  4.2.  In  this  case,  the  T-mesh  elements  are 
simply  subdivided.  Once  the  selected  T-mesh  elements  are  subdivided.  Steps  2  through  5  of  the 
local  refinement  algorithm  presented  in  Section  4.2  are  applied.  This  generates  the  final  refined 
analysis-suitable  T-spline  space.  Bezier  extraction  is  then  performed  resulting  in  a  new  set  of 
Bezier  elements.  This  process  is  then  repeated  until  the  resolution  of  the  T-spline  is  sufficient  for 
the  application. 

The  T-spline  container  ship  hull  in  Figure  18  is  first  designed  using  the  T-spline  plugin  for 
Rhino3d  [56],  It  is  then  transformed  into  an  analysis-suitable  model  using  local  refinement.  T- 
splines  are  a  popular  technology  in  ship  hull  design  because  an  entire  hull  can  be  modeled  by  a 
single  watertight  surface  with  a  minimal  number  of  control  points  [57].  In  analysis,  however,  far 
more  degrees-of-freedom  are  often  required  to  capture  the  physical  phenomenon  of  interest.  In 
other  words,  the  initial  T-mesh  must  undergo  additional  refinements  to  create  models  that  satisfy 
the  needs  of  analysis. 
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Figure  17:  An  analysis-suitable  local  refinement  framework. 
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Figure  18:  A  T-spline  container  ship  hull.  The  surface  is  C2 -continuous  everywhere. 


Figure  19:  The  regions  of  the  container  ship  hull  where  analysis-suitable  local  refinement  will  be  performed.  First, 
refinement  will  be  performed  in  the  rectangular  region  followed  by  highly  localized  refinement  along  the  curve. 
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To  demonstrate  the  pertinent  ideas,  we  assume  that  the  final  analysis-suitable  model  of  the 
hull  must  be  sufficiently  resolved  to  capture  the  response  of  the  ship  in  the  two  regions  outlined 
in  Figure  19.  This  will  require  refinements  in  the  region  of  the  hull  corresponding  to  the  rect¬ 
angle  followed  by  highly  localized  refinements  along  the  region  corresponding  to  the  curve.  Six 
iterations  of  refinement  are  performed  as  shown  in  Figures  20  through  25.  The  initial  T-spline  of 
the  hull  contains  just  75  control  points  and  36  Bezier  elements. 

The  first  iteration  of  local  refinement  of  the  ship  hull  is  shown  in  Figure  20.  A  set  of  Bezier 
elements  is  selected  for  refinement  as  shown  on  the  top  of  Figure  20.  The  refinement  framework 
described  in  Figure  17  is  then  applied.  First,  the  selected  Bezier  elements  are  used  to  select 
corresponding  T-mesh  elements.  These  T-mesh  elements  are  subdivided  and  the  local  refinement 
algorithm  described  in  Section  4.2  is  applied  to  the  resulting  subdivided  T-mesh.  The  control 
points  added  during  local  refinement  are  shown  in  the  middle  of  Figure  20.  Notice  that  these 
control  points  remain  localized  to  the  region  of  selected  Bezier  elements.  The  refined  control 
mesh  is  shown  on  the  bottom  of  Figure  20  with  the  new  control  points  highlighted.  The  refined 
set  of  Bezier  elements  is  then  extracted  from  the  refined  T-mesh  as  shown  in  Figure  21  on  the 
top.  We  note  that  the  transpose  of  the  refinement  operator,  MT  (see  Sections  2.4.2  and  4.3),  is 
used  to  update  control  point  positions  after  each  refinement  step.  This  ensures  that  the  geometry 
and  its  parameterization  are  preserved  exactly. 

The  next  five  iterations  of  local  refinement  are  shown  in  Figures  21  through  25.  In  Figure  21 
and  Figure  22  the  rectangular  region  in  Figure  19  undergoes  additional  local  refinement.  In 
Figures  23  through  25  highly  localized  refinement  is  performed  along  the  curve  in  Figure  19. 
Notice  that  the  refinement  pattern  follows  the  curve  without  excessive  propagation  of  control 
points  while  preserving  the  C~ -continuous  analysis-suitable  T-spline  basis.  The  fully  resolved 
analysis-suitable  model  has  1722  control  points  and  1922  Bezier  elements.  The  final  Bezier 
mesh  is  shown  in  Figure  25  on  the  bottom  and  the  final  T-mesh  is  shown  immediately  above 
it.  The  sequence  of  C2 -continuous  T-spline  spaces  is  nested  and  the  initial  geometry  is  exactly 
preserved  throughout. 

6.  Conclusion 

Analysis-suitable  T-splines  address  mathematical  and  practical  shortcomings  observed  when 
general  T-spline  spaces  are  used  as  a  basis  for  isogeometric  analysis.  Specifically,  analysis- 
suitable  T-splines  are  linearly  independent,  form  a  partition  of  unity,  and  can  be  locally  refined 
without  excessive  propagation.  We  have  developed  an  analysis-suitable  local  refinement  algo¬ 
rithm  for  this  class  of  T-splines,  and  implemented  it  for  bicubic  T-spline  surfaces. 

We  have  also  developed  an  efficient  adaptive  framework  which  combines  analysis-suitable 
T-splines,  Bezier  extraction,  and  analysis-suitable  local  refinement.  The  sequence  of  refined 
spaces  is  nested  and  exactly  preserves  the  initial  geometry.  The  basis  of  the  refined  T-spline 
spaces  maintains  the  smoothness  of  the  initial  basis.  We  have  demonstrated  its  effectiveness 
on  a  real-world  example  of  a  ship  hull  design.  The  procedures  described  provide  a  powerful 
methodology  for  instantiating  the  vision  of  isogeometric  analysis.  In  future  work,  we  plan  to 
describe  the  analysis-suitable  treatment  of  “extraordinary  points.” 
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Figure  20:  The  first  iteration  of  analysis-suitable  local  refinement  of  the  container  ship  hull  in  Figure  18.  The  Bezier 
elements  are  shown  on  the  top  with  those  elements  selected  for  refinement  highlighted  in  red.  The  control  points  added 
during  analysis-suitable  local  refinement  are  shown  in  the  middle.  The  refined  T-mesh  is  shown  on  the  bottom.  The  new 
control  points  are  highlighted  in  red.  (Note:  Some  of  the  T-mesh  edges  are  hidden  behind  the  ship  hull.) 
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Figure  21:  The  second  iteration  of  analysis-suitable  local  refinement  of  the  container  ship  hull  in  Figure  18.  The  Bezier 
elements  are  shown  on  the  top  with  those  elements  selected  for  refinement  highlighted  in  red.  The  control  points  added 
during  analysis-suitable  local  refinement  are  shown  in  the  middle.  The  refined  T-mesh  is  shown  on  the  bottom.  The  new 
control  points  are  highlighted  in  red.  (Note:  Some  of  the  T-mesh  edges  are  hidden  behind  the  ship  hull.) 
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Figure  22:  The  third  iteration  of  analysis-suitable  local  refinement  of  the  container  ship  hull  in  Figure  18.  The  Bezier 
elements  are  shown  on  the  top  with  those  elements  selected  for  refinement  highlighted  in  red.  The  control  points  added 
during  analysis-suitable  local  refinement  are  shown  in  the  middle.  The  refined  T-mesh  is  shown  on  the  bottom.  The  new 
control  points  are  highlighted  in  red.  (Note:  Some  of  the  T-mesh  edges  are  hidden  behind  the  ship  hull.) 
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Figure  23:  The  fourth  iteration  of  analysis-suitable  local  refinement  of  the  container  ship  hull  in  Figure  18.  The  Bezier 
elements  are  shown  on  the  top  with  those  elements  selected  for  refinement  highlighted  in  red.  The  control  points  added 
during  analysis-suitable  local  refinement  are  shown  in  the  middle.  The  refined  T-mesh  is  shown  on  the  bottom.  The 
new  control  points  are  highlighted  in  red.  Notice  that  analysis-suitable  local  refinement  remains  localized  to  the  Bezier 
elements  selected  for  refinement.  (Note:  Some  of  the  T-mesh  edges  are  hidden  behind  the  ship  hull.) 
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Figure  24:  The  fifth  iteration  of  analysis-suitable  local  refinement  of  the  container  ship  hull  in  Figure  18.  The  Bezier 
elements  are  shown  on  the  top  with  those  elements  selected  for  refinement  highlighted  in  red.  The  control  points  added 
during  analysis-suitable  local  refinement  are  shown  in  the  middle.  The  refined  T-mesh  is  shown  on  the  bottom.  The 
new  control  points  are  highlighted  in  red.  Notice  that  analysis-suitable  local  refinement  remains  localized  to  the  Bezier 
elements  selected  for  refinement.  (Note:  Some  of  the  T-mesh  edges  are  hidden  behind  the  ship  hull.) 
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Figure  25:  The  sixth  and  final  iteration  of  analysis- suitable  local  refinement  of  the  container  ship  hull  in  Figure  18.  The 
Bezier  elements  are  shown  on  the  top  with  those  elements  selected  for  refinement  highlighted  in  red.  The  control  points 
added  during  analysis-suitable  local  refinement  are  shown  below  that.  The  refined  final  T-mesh  is  then  shown.  The  new 
control  points  are  highlighted  in  red.  (Note:  Some  of  the  T-mesh  edges  are  hidden  behind  the  ship  hull.)  The  final  Bezier 
element  mesh  is  shown  on  the  bottom.  The  refinements  form  a  nested  sequence  of  C2  -continuous  spline  spaces.  The 
geometry  of  the  hull  is  unchanged  during  the  refinements. 


32 


Acknowledgements 

This  work  was  supported  by  grants  from  the  Office  of  Naval  Research  (N00014-08-C-0920 
and  N000 14-08- 1  -0992),  the  National  Science  Foundation  (CMI-0700807),  SINTEF  (UTA10- 
000374),  the  National  Science  Foundation  of  China  (11031007  and  60903148),  the  Chinese 
Universities  Scientific  Fund,  and  the  Chinese  Academy  of  Science  (Startup  Scientific  Research 
Foundation).  M.A.  Scott  was  partially  supported  by  an  ICES  CAM  Graduate  Fellowship.  This 
support  is  gratefully  acknowledged. 


References 

[1]  T.  J.  R.  Hughes,  J.  A.  Cottrell,  Y.  Bazilevs,  Isogeometric  analysis:  CAD,  finite  elements,  NURBS,  exact  geometry, 
and  mesh  refinement,  Computer  Methods  in  Applied  Mechanics  and  Engineering  194  (2005)  4135 — 4195. 

[2]  J.  A.  Cottrell,  T.  J.  R.  Hughes,  Y.  Bazilevs,  Isogeometric  analysis:  Toward  Integration  of  CAD  and  FEA,  Wiley, 
Chichester,  2009. 

[3]  J.  A.  Cottrell,  A.  Reali,  Y.  Bazilevs,  T.  J.  R.  Hughes,  Isogeometric  analysis  of  structural  vibrations,  Computer 
Methods  in  Applied  Mechanics  and  Engineering  195  (2006)  5257-5296. 

[4]  J.  A.  Evans,  Y.  Bazilevs,  I.  Babuska,  T.  J.  R.  Hughes,  rc-widths,  sup-infs,  and  optimality  ratios  for  the  ^-version  of 
the  isogeometric  finite  element  method,  Computer  Methods  in  Applied  Mechanics  and  Engineering  198  (21-26) 
(2009)  1726-1741. 

[5]  I.  Akkerman,  Y.  Bazilevs,  V.  Calo,  T.  J.  R.  Hughes,  S.  Hulshoff,  The  role  of  continuity  in  residual-based  variational 
multiscale  modeling  of  turbulence.  Computational  Mechanics  41  (2008)  371-378. 

[6]  Y.  Bazilevs,  V.  M.  Calo,  J.  A.  Cottrell,  T.  J.  R.  Hughes,  A.  Reali,  G.  Scovazzi,  Variational  multiscale  residual-based 
turbulence  modeling  for  large  eddy  simulation  of  incompressible  flows,  Computer  Methods  in  Applied  Mechanics 
and  Engineering  197  (2007)  173-201. 

[7]  Y.  Bazilevs,  I.  Akkerman,  Large  eddy  simulation  of  turbulent  Taylor-Couette  flow  using  isogeometric  analysis  and 
residual-based  variational  multiscale  method,  Journal  of  Computational  Physics  229  (2010)  3402-3414. 

[8]  Y.  Bazilevs,  C.  Michler,  V.  M.  Calo,  T.  J.  R.  Hughes,  Isogeometric  variational  multiscale  modeling  of  wall-bounded 
turbulent  flows  with  weakly  enforced  boundary  conditions  on  unstretched  meshes,  Computer  Methods  in  Applied 
Mechanics  and  Engineering  199  (13-16)  (2010)  780-790. 

[9]  Y.  Bazilevs,  V.  M.  Calo,  Y.  Zhang,  T.  J.  R.  Hughes,  Isogeometric  fluid-structure  interaction  analysis  with  applica¬ 
tions  to  arterial  blood  flow,  Computational  Mechanics  38  (2006)  310-322. 

[10]  Y.  Bazilevs,  V.  M.  Calo,  T.  J.  R.  Hughes,  Y.  Zhang,  Isogeometric  fluid- structure  interaction:  Theory,  algorithms, 
and  computations,  Computational  Mechanics  43  (2008)  3-37. 

[11]  Y.  Zhang,  Y.  Bazilevs,  S.  Goswami,  C.  Bajaj,  T.  J.  R.  Hughes,  Patient-specific  vascular  NURBS  modeling  for 
isogeometric  analysis  of  blood  flow,  Computer  Methods  in  Applied  Mechanics  and  Engineering  196  (2007)  2943- 
2959. 

[12]  Y.  Bazilevs,  J.  R.  Gohean,  T.  J.  R.  Hughes,  R.  D.  Moser,  Y.  Zhang,  Patient-specific  isogeometric  fluid-structure 
interaction  analysis  of  thoracic  aortic  blood  flow  due  to  implantation  of  the  Jarvik  2000  left  ventricular  assist 
device.  Computer  Methods  in  Applied  Mechanics  and  Engineering  198  (45-46)  (2009)  3534-3550. 

[13]  F.  Auricchio,  L.  B.  da  Veiga,  C.  Lovadina,  A.  Reali,  The  importance  of  the  exact  satisfaction  of  the  incompress¬ 
ibility  constraint  in  nonlinear  elasticity:  mixed  FEMs  versus  NURBS-based  approximations,  Computer  Methods 
in  Applied  Mechanics  and  Engineering  199  (2010)  314-323. 

[14]  F.  Auricchio,  L.  B.  da  Veiga,  A.  Buffa,  C.  Lovadina,  A.  Reali,  G.  Sangalli,  A  fully  “locking-free”  isogeomet¬ 
ric  approach  for  plane  linear  elasticity  problems:  A  stream  function  formulation,  Computer  Methods  in  Applied 
Mechanics  and  Engineering  197  (2007)  160-172. 

[15]  T.  Elguedj,  Y.  Bazilevs,  V.  M.  Calo,  T.  J.  R.  Hughes,  B  and  F  projection  methods  for  nearly  incompressible 
linear  and  non-linear  elasticity  and  plasticity  using  higher-order  NURBS  elements,  Computer  Methods  in  Applied 
Mechanics  and  Engineering  197  (2008)  2732-2762. 

[16]  J.  A.  Cottrell,  T.  J.  R.  Hughes,  A.  Reali,  Studies  of  refinement  and  continuity  in  isogeometric  analysis.  Computer 
Methods  in  Applied  Mechanics  and  Engineering  196  (2007)  4160^-183. 

[17]  D.  J.  Benson,  Y.  Bazilevs,  M.  C.  Hsu,  T.  J.  R.  Hughes,  Isogeometric  shell  analysis:  The  Reissner-Mindlin  shell, 
Computer  Methods  in  Applied  Mechanics  and  Engineering  199  (5-8)  (2010)  276  -  289. 

[18]  D.  J.  Benson,  Y.  Bazilevs,  M.  C.  Hsu,  T.  J.  R.  Hughes,  A  large  deformation,  rotation-free,  isogeometric  shell, 
International  Journal  for  Numerical  Methods  in  Engineering,  in  press,  doi:  10. 10 16/j.cma.2010. 12.003. 


33 


[19]  D.  J.  Benson,  Y.  Bazilevs,  E.  De  Luycker,  M.  C.  Hsu,  M.  A.  Scott,  T.  J.  R.  Hughes,  T.  Belytschko,  A  generalized 
finite  element  formulation  for  arbitrary  basis  functions:  from  isogeometric  analysis  to  XFEM,  International  Journal 
for  Numerical  Methods  in  Engineering  83  (2010)  765-785. 

[20]  R.  Echter,  M.  Bischoff,  Numerical  efficiency,  locking  and  unlocking  of  NURBS  finite  elements,  Computer  Methods 
in  Applied  Mechanics  and  Engineering  199  (5-8)  (2010)  374-382. 

[21]  J.  Kiendl,  Y.  Bazilevs,  M.  C.  Hsu,  R.  Wuechner,  K.  U.  Bletzinger,  The  bending  strip  method  for  isogeometric  anal¬ 
ysis  of  Kirchhoff-Love  shell  structures  comprised  of  multiple  patches,  Computer  Methods  in  Applied  Mechanics 
and  Engineering  199  (37-40)  (2010)  2403-2416. 

[22]  H.  Gomez,  V.  M.  Calo,  Y.  Bazilevs,  T.  J.  R.  Hughes,  Isogeometric  analysis  of  the  Cahn-Hilliard  phase-field  model, 
Computer  Methods  in  Applied  Mechanics  and  Engineering  197  (2008)  4333^-352. 

[23]  H.  Gomez,  T.  J.  R.  Hughes,  X.  Nogueira,  V.  M.  Calo,  Isogeometric  analysis  of  the  isothermal  Navier-Stokes- 
Korteweg  equations.  Computer  Methods  in  Applied  Mechanics  and  Engineering  199  (25-28)  (2010)  1828-1840. 

[24]  S.  Lipton,  J.  A.  Evans,  Y.  Bazilevs,  T.  Elguedj,  T.  J.  R.  Hughes,  Robustness  of  isogeometric  structural  discretizations 
under  severe  mesh  distortion,  Computer  Methods  in  Applied  Mechanics  and  Engineering  199  (5-8)  (2010)  357  - 
373. 

[25]  W.  A.  Wall,  M.  A.  Frenzel,  C.  Cyron,  Isogeometric  structural  shape  optimization,  Computer  Methods  in  Applied 
Mechanics  and  Engineering  197  (2008)  2976-2988. 

[26]  X.  Qian,  Full  analytical  sensitivities  in  NURBS  based  isogeometric  shape  optimization,  Computer  Methods  in 
Applied  Mechanics  and  Engineering  199  (29-32)  (2010)  2059-2071. 

[27]  A.  P.  Nagy,  M.  M.  Abdalla,  Z.  Gurdal,  Isogeometric  sizing  and  shape  optimisation  of  beam  structures.  Computer 
Methods  in  Applied  Mechanics  and  Engineering  199  (17-20)  (2010)  1216-1230. 

[28]  A.  P.  Nagy,  M.  M.  Abdalla,  Z.  Gurdal,  On  the  variational  formulation  of  stress  constraints  in  isogeometric  design, 
Computer  Methods  in  Applied  Mechanics  and  Engineering  199  (41-44)  (2010)  2687-2696. 

[29]  A.  Buffa,  G.  Sangalli,  R.  Vazquez,  Isogeometric  analysis  in  electromagnetics:  B-splines  approximation,  Computer 
Methods  in  Applied  Mechanics  and  Engineering  199  (17-20)  (2010)  1143-1152. 

[30]  E.  Cohen,  T.  Martin,  R.  M.  Kirby,  T.  Lyche,  R.  F.  Riesenfeld,  Analysis-aware  modeling:  Understanding  quality 
considerations  in  modeling  for  isogeometric  analysis,  Computer  Methods  in  Applied  Mechanics  and  Engineering 
199  (5-8)  (2010)  334—356. 

[31]  J.  Lu,  Circular  element:  Isogeometric  elements  of  smooth  boundary,  Computer  Methods  in  Applied  Mechanics  and 
Engineering  198  (30-32)  (2009)  2391-2402. 

[32]  H.-J.  Kim,  Y.-D.  Seo,  S.-K.  Youn,  Isogeometric  analysis  for  trimmed  CAD  surfaces,  Computer  Methods  in  Applied 
Mechanics  and  Engineering  198  (37-40)  (2009)  2982-2995. 

[33]  P.  Costantini,  C.  Manni,  F.  Pelosi,  M.  L.  Sampoli,  Quasi-interpolation  in  isogeometric  analysis  based  on  generalized 
B-splines,  Computer  Aided  Geometric  Design  27  (8)  (2010)  656-668. 

[34]  M.  Aigner,  C.  Heinrich,  B.  Juettler,  E.  Pilgerstorfer,  B.  Simeon,  A.  V.  Vuong,  Swept  Volume  Parameterization  for 
Isogeometric  Analysis,  in:  Hancock,  E.  R.  and  Martin,  R.  R.  and  Sabin,  M.  A.  (Ed.),  Mathematics  of  Surfaces  XIII, 
Vol.  5654  of  Lecture  Notes  in  Computer  Science,  2009,  pp.  19-44. 

[35]  T.  Martin,  E.  Cohen,  R.  M.  Kirby,  Volumetric  parameterization  and  trivariate  B-spline  fitting  using  harmonic  func¬ 
tions,  Computer  Aided  Geometric  Design  26  (6)  (2009)  648-664. 

[36]  W.  Wang,  Y.  Zhang,  Wavelets-based  NURBS  simplification  and  fairing,  Computer  Methods  in  Applied  Mechanics 
and  Engineering  199  (5-8)  (2010)  290-300.  doi:10.1016/j.cma.2009.04.003. 

[37]  W.  Wang,  Y.  Zhang,  M.  A.  Scott,  T.  J.  R.  Hughes,  Converting  an  unstructured  quadrilateral  mesh  to  a  standard 
T-spline  surface,  Computational  Mechanics,  accepted  for  publication. 

[38]  T.  W.  Sederberg,  J.  Zheng,  A.  Bakenov,  A.  Nasri,  T-splines  and  T-NURCCSs,  ACM  Transactions  on  Graphics  22 
(3) (2003)  477-484. 

[39]  T.  W.  Sederberg,  G.  T.  Finnigan,  X.  Li,  H.  Lin,  Watertight  trimmed  NURBS,  ACM  Transactions  on  Graphics  27  (3) 
(2008)  Article  no.  79. 

[40]  H.  Ipson,  T-spline  merging.  Master’s  thesis,  Brigham  Young  University  (April  2005). 

[41]  T.  W.  Sederberg,  D.  L.  Cardon,  G.  T.  Finnigan,  N.  S.  North,  J.  Zheng,  T.  Lyche,  T-spline  simplification  and  local 
refinement,  ACM  Transactions  on  Graphics  23  (3)  (2004)  276-283. 

[42]  Y.  Bazilevs,  V.  M.  Calo,  J.  A.  Cottrell,  J.  A.  Evans,  T.  J.  R.  Hughes,  S.  Lipton,  M.  A.  Scott,  T.  W.  Sederberg, 
Isogeometric  analysis  using  T-splines,  Computer  Methods  in  Applied  Mechanics  and  Engineering  199  (5-8)  (2010) 
229-263. 

[43]  M.  Dorfel,  B.  Jiittler,  B.  Simeon,  Adaptive  isogeometric  analysis  by  local  /z-refinement  with  T-splines,  Computer 
Methods  in  Applied  Mechanics  and  Engineering  199  (5-8)  (2009)  264-275. 

[44]  M.  J.  Borden,  M.  A.  Scott,  J.  A.  Evans,  T.  J.  R.  Hughes,  Isogeometric  finite  element  data  structures  based 
on  Bezier  extraction  of  NURBS,  International  Journal  for  Numerical  Methods  in  Engineering,  in  press,  doi: 
10.1002/nme.2968. 

[45]  M.  A.  Scott,  M.  J.  Borden,  C.  V.  Verhoosel,  T.  J.  R.  Hughes,  Isogeometric  Finite  Element  Data  Structures  based 

34 


on  Bezier  Extraction  of  T-splines,  International  Journal  for  Numerical  Methods  in  Engineering,  in  press,  doi: 
10.1 002/nme.3 1 67. 

[46]  C.  V.  Verhoosel,  M.  A.  Scott,  R.  de  Borst,  T.  J.  R.  Hughes,  An  isogeometric  approach  to  cohesive  zone  modeling, 
International  Journal  for  Numerical  Methods  in  Engineering,  in  press,  doi:  10.1002/nme.3061. 

[47]  C.  V.  Verhoosel,  M.  A.  Scott,  T.  J.  R.  Hughes,  de  Borst,  R.,  An  isogeometric  analysis  approach  to  gradient  damage 
models.  International  Journal  for  Numerical  Methods  in  Engineering,  in  press,  doi:  10.1002/nme.3150. 

[48]  X.  Li,  J.  Zheng,  T.  W.  Sederberg,  T.  J.  R.  Hughes,  M.  A.  Scott,  On  the  linear  independence  of  T- 
splines.  Computer  Aided  Geometric  Design,  submitted  for  publication  (see  also  ICES  Report  10-40  at 
http :  //www.  ice  s .  utexas .  edu/re  search/reports ) . 

[49]  X.  Li,  J.  Zheng,  T.  W.  Sederberg,  On  T-spline  classification,  in  preparation. 

[50]  G.  T.  Finnigan,  Arbitrary  degree  T-splines,  Master’s  thesis,  Brigham  Young  University  (August  2008). 

[51]  T.  W.  Sederberg,  J.  Zheng,  X.  Song,  Knot  intervals  and  multi-degree  splines.  Computer  Aided  Geometric  Design 
20  (2003)  455^168. 

[52]  T.  J.  R.  Hughes,  The  Finite  Element  Method:  Linear  Static  and  Dynamic  Finite  Element  Analysis,  Dover  Publica¬ 
tions,  Mineola,  NY,  2000. 

[53]  A.  Buffa,  D.  Cho,  G.  Sangalli,  Linear  independence  of  the  T-spline  blending  functions  associated  with  some  par¬ 
ticular  T-meshes,  Computer  Methods  in  Applied  Mechanics  and  Engineering  199  (23-24)  (2010)  1437-1445. 

[54]  X.  Li,  M.  A.  Scott,  The  dimension  of  spline  spaces  over  the  extended  T-mesh,  in  preparation. 

[55]  M.  J.  Borden,  M.  A.  Scott,  C.  V.  Verhoosel,  C.  M.  Landis,  T.  J.  R.  Hughes,  Phase-field  modeling  of  dynamic 
fracture  using  isogeometric  analysis,  in  preparation. 

[56]  T-Splines,  Inc.,  http://www.tsplines.com/rhino/  (2011). 

[57]  M.  T.  Sederberg,  T.  W.  Sederberg,  T-splines:  A  technology  for  marine  design  with  minimal  control  points,  in:  Pro¬ 
ceedings  of  the  Chesapeake  Powerboat  Symposium,  http://www.tsplines.com/rn/T-SplinesChesapeakepaper.pdf, 
2010. 


35 


